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Abstract. This work aims to extend part of the two dimensional results of 
Duplantier and Sheffield on Liouville quantum gravity |DS11) to four dimen- 
sions, and indicate possible extensions to other even-dimensional spaces R-^" 
as well as Riemannian manifolds. 

Let © be the Gaussian free field on R* with the underlying Hilbert space 
H-^ (R^) and the inner product ^(/ — A)^-,-^ ^, and 6 a generic element 

from 0. We consider a sequence of random Borel measures on R*, denoted 
by {mf^ (dx) : n > l}, each of which is absolutely continuous with respect to 
the Lebesgue measure dx, and the density function is given by the exponential 
of a centered Gaussian family parametrized by x £ R*. We show that with 
probability 1, mf^ (dx) weakly converges as en J, 0, and the limit measure can 
be "formally" written as "m" (dx) = e'^~>^^^^ dx". In this setting, we also prove 
a KPZ relation, which is the quadratic relation between the scaling exponent 
of a bounded Borel set on R** under the Lebesgue measure and its counterpart 
under the random measure (dx). 

Our approach is similar to the one used in [DSll] with adaptations to W^. 



1. Introduction 

Random measures have long been considered in 2-dimensional conformal field 
theory and quantum gravity since the work of Knizhnik, Polyakov and Zamolod- 
chikov |KPZ| . Recently, a probabilistic proof of the formula due to Knizhnik, 
Polyakov and Zamolodchikov was given by Duplantier and Shefheld in |DS11| . On 
the unit planar disc D, Duplantier and Sheffield construct the Liouville quantum 
gravity measure "e'^'^'^^-'dz", where dz is the Lebesgue measure on D, 7 is a properly 
chosen positive constant and h is an instance of the Gaussian free field (GFF) on 
D with the Dirichlet inner product. To be specific, they prove that the random 
measure exists as the weak convergence limit of /^e'^'^'^^^dz as e J, 0, where he {z) 
is the circular average of h over the circle centered at z with radius e. They further 
show that there is a quadratic relation, known as the KPZ relation, between the 
scaling exponent of a random set under the Lebesgue measure, and its counterpart 
under the quantum gravity measure. Another derivation was given in |DB| . Gauss- 
ian free field in dimension 2 has also been considered in |HMP| and numerous other 
papers. 

In this article, we generalize part of the results from |DS11| to four dimensions. 
We define the Euclidean GFF on K^, denoted by O, with the inner product deter- 
mined by the Bessel operator (/ — A)^. In other words, the underlying Hilbert space 
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of 6 is given by the Sobolev space (M'') with the inner product \ {I — A)^ •, • j 

In this setting, we prove (Section 2, Theorem [5]) that given < 7^ < 27r^, almost 
every 6* £ 8 admits a random measure on R** which "formally" has the density 
g27S(a;) with respect to the Lebesgue measure dx on M"*. We also show that this 
random measure satisfies a KPZ relation similar to the one in the two dimensional 
case. Namely, if k G [0, 1] is the scaling exponent of a bounded Borel set in 
under the Lebesgue measure, and K e [0, 1] is the scaling exponent of the same 
set but under the random measure obtained above (both k and K will be defined 
in Section 4), then k, and K satisfy the following quadratic relation (Section 4, 
Theorem HI): 

Our proof follows the outline of the proof in |DS11| with adaptations to four 
dimensions. Mainly we have to overcome (both in "designing" the model to work 
with and in technical details) the difficulties caused by the absence in our problem 
the two dimensional conformal structure. To interpret rigorously an instance 9 
of the GFF on the entire Euclidean space R^, we adopt the theory of the abstract 
Wiener space. A key ingredient in this theory is the underlying Hilbert space whose 
inner product determines the covariance structure of the field. It is already known 
that in order to obtain a measure which "formally" has the exponential of 9 [x) 
as the density with respect to dx, the covariance function Gov {9 (x) , 9 (y)) can 
at most grow at the rate of — log \x — y\ when \x — y\ is small. Taking this into 

account, iJ^ (R'') with the inner product |^(/ — A)^-,-^ ^ becomes our natural 

choice. Also this way of defining the GFF makes it possible, in certain situations, 
to obtain explicit formulas of the covariance function. To construct the random 
measure and thereafter to study it, we always need to relate it to a sequence of 
approximating measures which converges in some proper sense. So it is our intention 
to choose the approximating measures appropriately so they will be convenient to 
work with. In the two dimensional case in |DS11| . the approximating measures are 
in terms of the circular averages of the GFF on D. In fact, the properties of the 
Gaussian family consisting of these circular averages play an important role in the 
proof. For example, if h is an instance of the GFF on D, then given any z e D, 
the one-parameter family {h^ (z) : < e < 1} has up to a time change the same 
distribution as a standard Brownian motion. Such properties are derived from the 
Green's function of the Laplace operator A on D, which, in particular, is harmonic. 
Therefore, it should not be surprising that the trivial analogue in four dimensions, 
that is, the family of spherical averages of 9, fails to have such properties, which 
makes it a less than optimal substitute for in carrying out this project on R^. 
In Section 2, we present one possible replacement for he in four dimensions which 
still has simple and concrete geometric interpretations (in fact, it is given by a 
functional of the spherical average of 9), but possesses, to a large extent, similar 
properties to those of in two dimensions. In Section 3, we use the results from 
Section 2 to build the approximating measures, and then prove they almost surely 
admit a limit measure in the sense of weak convergence. In Section 4, we lay out an 
outline to derive the KPZ relation and the proofs of the main results are collected 
in Section 5. Other work on the KPZ relation in higher dimensions with different 
settings can be found in |JJRV| IRV] . 
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Our original interest in constructing such a random measure lies in its poten- 
tial applications in the study of conformal classes of Riemannian metrics. In fact, 
another more geometric point of view on the GFF on a planar domain or more 
generally on a surface S, is to consider it as a measure on a conformal class of 
metrics on E, where the measure is constructed with the help of a reference metric 
go on S, but where the result does not depend on go- It seems natural to generalize 
this approach to conformal classes of metrics on higher-dimensional manifolds. It 
turns out that on a compact four-dimensional manifold 971, a natural replacement 
for the Laplace-Beltrami operator A (that is used in the construction of the GFF on 
surfaces) is the 4-th order Paneitz operator that arises in the conformal geometry. 
More generally, on compact 2ri-dimensional manifolds it seems natural to use the 
dimension- critical GJMS operator in the construction of higher-dimensional ana- 
logues of the two-dimensional GFF. This will be further explained in the second 
part of Section 6. 

Acknowledgments. The authors thank L. Addario-Berry and D. Stroock for their 
valuable comments on an earlier version of this article. We also thank R. Adler, I. 
Binder, Y. Canzani, B. Duplantier, S. Klevtsov, R. Ponge, P. Sarnak, S. Sheffield, 
J. Taylor, I. Wigman, P. Yang, S. Zelditch for many fruitful conversations. 



We start with a brief review of some fundamental facts about the abstract Wiener 
space theory ( |Gr| or [S08J). An abstract Wiener space is commonly used in con- 
structions of infinite dimensional Gaussian measures. The basic setting of an ab- 
stract Wiener space is as follows. Given a (infinite dimensional) Banach space Q 
and a (infinite dimensional) Hilbert space H, assume both Q and H are separable, 
and H can be continuously embedded into as a dense subspace. Therefore if x* 
is a bounded linear functional on Q (denoted by x* € 6*), then there is unique 
hx* £ H such that for every h G H, {h, hx')^ — {h,x*), where (•, *) refers to the 
action of 9* on Q (or more specifically in later discussions, the action of tempered 
distributions on test functions). Let W be a probability measure on (0, Se) where 
Se is the Borel cr— algebra of 8. If W satisfies 



then the triple {H,Q,yV) is called an abstract Wiener space. It is known ([SllJ, 
§8.3) that given any separable Hilbert space, one can always find Q and W such 
that {H,Q,yV) forms an abstract Wiener space. Moreover, since {h^* : x* £ &*} 
is also dense in H, the linear mapping 



can be uniquely extended as a linear isometry from H to i^(yV). Its images 
{1 (h) : h £ H}, known as the Paley-Wiener integrals, form a centered Gaussian 
family whose covariance is given by 



We point out that although the Hilbert structure of H plays an essential role, 
W {H) = and the choice of Q is not unique. 



2. Spherical Averages of GFF on M"^ 




I: £ H ^I{hx') = {■,x*) £ (W) 



[I{hi)I{h2)] = {hi,h2)H for aU hi,h2 £ H. 
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As we have mentioned in the previous section, we consider in our project the in- 
finite dimensional Gaussian measure on the space of certain tempered distributions 
on M"*, with the underlying Hilbert space given by the Sobolev space H = (K**), 
which is the completion of the real valued Schwartz test function space S (M*) 
under the inner product 

(/i, f2)H = f {I- A)' /i (x) /2 (x) dx for aU /i, ^ G 5 (R^) . 

Then, given this particular choice of H, our notion of the Gaussian free field on 
refers to any probability space (8,*Be,VV) such that {Q,H,yV) forms an abstract 
Wiener space. For example, if Q is the space of continuous functions 9 : M* R 
satisfying 

lim (log(e + |6'(a;)| = 0, 

i 

then 9 can be chosen as the image of 6 under the Bessel operator (/ — A) * , i.e., 

e = |(/- A)3 6': e/e e}. 

From this we observe that Q consists of tempered distributions which in general 
are not defined point-wise. Nonetheless, we can understand some properties of the 
GFF by studying the Paley- Wiener integrals, which can be viewed as "generalized" 
action of certain tempered distributions on Q. 

In addition, if H^^ = H^^ (R'') is the Hilbert space consisting of tempered 
distributions /i such that 

IImII^-. ^ / (i + |er)"'lA(Ol'rfe<oo 

where /t is the Fourier transform (without the factor (27r) ^ in the definition) of /i, 
then we can identify H with since (/ — A) : H^^ — > is obviously a linear 
isometry. We will abuse the notatiorQby denoting "/ii/" the image of G H^^ under 
(/- A)" . Then is the unique element in H such that {h, v) — {h, hif)jj for all 
h ^ H, which suggests that the corresponding Paley- Wiener integral I{hif) can be 
viewed as a "representation" of the action of on 0, even though v is not in 8* and 
X{h^) {9) is only defined for almost every e 8. Meanwhile, : u E iJ^^} 

is also a Gaussian family whose covariance is given by 

With these in mind, as a natural analogue of the 2D circular average, we consider 
the spherical average of the GFF on R''. To this end, for every a; e R'' and e > 0, 
denote the tempered distribution determined by 

(/, a^) s [ f (y) da (y) for aU / G 5 (R^) , 

where Se (x) is the sphere centered at x with radius e, and da is the surface area 
measure on Se (x). Clearly, the action of cr^ is to take the spherical average of / 



^The subscript of "/ij/" is an element of H ^, not to be confused with "h^;*" in the definition 
of the abstract Wiener space where x* € 0*. 
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over Se (x). It is an easy matter to verify that a'f ^ H ^. In fact, one only needs 
to write down the Fourier transform of cr^ as 

(2.1) al(0=2(6|C|)-'Ji(e|^|)e^(-'«V 

where Jk (r) is the Bessel function of order fc e N, and use the fact that Jk (r) is 
asymptotic to r~^/^ when r is large. 

As we have indicated in the introduction (and as we will confirm in the next 
lemma), the spherical average of the GFF on does not behave as "nicely" as the 
circular average of the GFF in two dimensions. For one thing : e > O} 

fails to be a reversed Markov process. An intuitive way to view this is that, the 
spherical average does not bear enough information in itself for this Gaussian fam- 
ily parametrized by radius e > to be (reversed) Markovian. It might also be 
helpful to relate this to the following analogous problem: when solving PDEs with 
higher order differential operator on a domain with boundary, one often needs more 
than one boundary condition (e.g., both the Dirichlet and the Neumann boundary 
conditions) to uniquely determine the solution. Inspired by this idea, besides the 
average itself we will also "collect" one more piece of information about the GFF 
from each sphere, which is the "derivative" of the average with respect to the radius. 
Namely, for every a; G M'' and e > 0, denote da^ the tempered distribution given 
by (/, dcTg ) = ^ (/, (T^ ) for all / G 5 , then the action of da^ can be viewed as 
to take the derivative of the spherical average of the GFF in the radial direction. 
It follows trivially from (|2.ip that 

(2.2) daf (C) = (0 - -2e-'j2 (e 1^1) e^(^'«)«^ 

In particular, da^ is also in and so {l (ft-ay ) ,1 (hda^) : a; € R'*, e > O} forms a 
centered Gaussian family whose covariance is determined by the inner product 
of {a^,da'^ : cc e R'i,e > O}. 

The next lemma in some sense validates our decision to take da^ into account. 
It shows that by putting I (ha^) and its "derivative" X (hda^) togethei0, not only 
does the Gaussian family recover the reversed Markov property in the concentric 
case (with x € fixed, parametrized by e > only), the non-concentric family 
(parametrized by both e > and x S M"*) also resembles, to a large extent, its 
counterpart in two dimensions. To be precise, we define the vector-valued Gaussian 
random variable: 

= ( ^il^^'l ^ for every x £ and e > 0. 

Then, under certain circumstances, the covariance matrix of the Gaussian family 
can be evaluated explicitly as follows. 

Lemma 1. For r e (0, oo), define the following four matrices: 

"^^"^'-y K'{ (r) -K2 (r) /r )''^'^^>-\h (r) /r /(' (r) 



h [r) h{r)/r r ^ ' ' \ [r) /r 



^This idea came from discussions with Daniel W. Stroock when the first author was studying 



at MIT. 
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where Ik,Kk are the modified Bess el functions of order fc € N. Then, 
(1), given a; G and ei > £2 > 0, 



(2.3) 



E 



47r2) 



A(ei)B^(e2) 



In particular, {V^ : e > 0} is a vector-valued Gaussian reversed Markov process in 
the sense that for every Borel A C R^, 

where a \ : ?7 > ei} is the a— algebra generated by {V^ : > ci}- 
(2), given x,y £ , x ^ y, and ei, e2 > with ei> \x — y\ + £2; 



(2.4) 



E 



27r' 



2 , A(ei)C(|x-y|)B' (£2) 



(3), given x,y G M.'^, x ^ y, and £1, £2 > with \x — y\ > £1 + £2, 



(2.5) 



E 



w 



T 



27r 



2 ,B(£i)D(|x-2/|)BT(£2). 



The proof of (|2.3p - (|2.5p relies heavily on the integral formulas and identities of 
Bessel functions. The complete detailed computations are given in the appendix. 
Here we only make the following observations. 

First, we claim that the distribution of {V/ : a; G M^,£ > 0} is invariant under 

isometries in spatial variables in the sense that {y7("^xGR4,£>0} has exactly 

the same distribution as {V/ : a; G M'',£ > O} for every T : M** — M"* satisfying 
\T {x) — T(j/)| = — 2/1 for all x,y E M**. Perhaps the most straightforward way 
to see this is to write down the covariance matrix of the family, or equivalently, 
the inner product of {af, daf : x G M"*, £ > O} in the integral form using ()2.ip 
and p.2p . We will actually do this in the appendix (formulas (|7.ip - (|7.6p ). The 
result shows that the only dependence of the covariance matrix on spatial variables 
is through the distance between centers of the spheres that are involved. 

Second, we point out that all the matrices above A (r), B (r), C (r) and D (r) are 
invertible for all r > 0. This fact can certainly be verified by direct computations 
using the explicit formulas given above, but it also follows, more generally, from 
the simple fact that da^ is linearly independent of for every x,y G and 



£ > 0. Therefore, assuming (|2.3p is true, then given x fixed and £1 > £2 > 0, the 
conditional expectation of conditioning on V^^ equals 



E 



w 



yx 



B(£2)B-1(£i)K:- 



On the other hand, we observe that for all ?y > £1, 



E 



w 



This means, 



(K:-B(£2)B-l(£l)K:)ra" 
K^-B(£2)B-1(£i)F,^ 



= 0. 



is independent of V^^ for all 77 > £1, which certainly implies the reversed Markov 
property. 
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Next, we observe that under the circumstances as prescribed in Lemma [1] the 
covariance matrix of {Vf : a; e R"*, e > O} is "separable" in the sense that it sphts 
into factors each of which only depends on one of the variables ei, 62 and — yj. 
A second look at the formulas (|2.3p - (l2.5p suggests that we should "normalize" 
by (e). Namely, if denote [/f = (e) V/, then the previous observations 
imply that given x fixed, {U^ : e > 0} is a vector-valued Gaussian process with 
independent (reversed) increment whose distribution does not depend on x. More- 



over, 



and show that E 



w 



only depends on ei and |a; — 



whei|^ -Bej (y) ^ (^)i Btiid the same matrix only depends on \x — y\ when 
{y) n -Bfi (x) = 0. Because C/f has these properties, we are one step closer 
to finding a plausible replacement for the circular average of the two dimensional 
GFF. 

Clearly, for any constant C = (Ci, C2) € R^^ {U^ X)m^ 

is a scalar valued Gaussian 

random variable (in fact, it is a Paley- Wiener integral), which, when parametrized 
by a; € and e > 0, forms a Gaussian family that preserves the properties de- 
scribed above. Our goal is to find a proper C G such that the random variable 

becomes a "legitimate" approximation for a multiple of the value of 9 at point x for 
every a; S M^. Namely, we want to choose C so that if /ij S H^^ is given by 



(2.6) 



(e) 



at 



then /ig converges to a constant multiple of the point mass (^j, at x as e 4, in the 
sense of tempered distribution. We can reach this goal by writing down the formula 
of B^^ (e) explicitly and examining the asymptotics of the Bessel functions near 
the origin (detailed computations are given in the appendix). As a result, we find 
that C = (1, 1)^ will serve the purpose, in which case /ij — > 2Sx as e J, for every 
x G R'*. From now on, we will assume /x^ is as in p.6p with C — (1, 1)^. Since 
= ^^B~^ (e) V"/, we can transfer the results in Lemma[T]to the Gaussian 
family {I (/i^j) : x G M^, e > O}. 

Theorem 2. Define the positive function G : r G (0, 00) 1— > G (r) G (0, 00) by 

1 \ 2/1 (r) Ki (r) + 2/2 (r) Ko (r) - 1 



47r2 



(2.7) G(r)E 
Then, we have 

(1), given a; G M'* and ei > £2 > 0, 



Il{r)-I,{r)h{r) 



(2.8) 



lih 



I h 



E 



w 



(/i^.J] -G(ei), 



In particular, \l (^m?) ^ e > O} is a Gaussian process with independent reversed in- 
crements in the sense thati {hfi^ ^ —I (j^fj.^ ) '■^ independent of a |z (j^fi^j : 77 > ei|. 



(2), given y G 



, X ^ y, and ei, e2 > with ei > |a; — y\ 

1 



(2.9) E^ [l(/i^.Jl(/i^«J] ^h{\x-y\)G{e^) 



£2, 
/2 



4^2/2(£^)„/Q(,^)/2(e^)- 



■^The notation "B^ (i^)" (i^)") denotes the open (closed) ball centered at x with radius 
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(3), given x,y £ , x ^ y, and ei, 62 > with |x — 2;| > ei + £2, 

(2.10) w}^\x(h^Ax(h^A] = -Li^o (|x - yl) . 



There is not much to be said about the proof since everything follows from 
straightforward computations and the results in Lemma [TJ However, we should 
point out the following facts. 

First, it is an easy matter to check that G : (0,oo) (0,oo) is smooth and 
strictly decreasing with limr^o G (r) = +00 and limrfoo G {r) — 0. Therefore, G~^ 
is defined and also strictly decreasing on (0, 00). Fix a positive constant i?, and for 
every r G (0, i?], define 

(2.11) 0<tEEG(r)-G(i?) andXt=l(v , 

Then {Xt : t > 0}, as a Gaussian process on 8 under W, has the same distribution 
as the standard Brownian motion, which in particular is independent of the choice 
of X. In other words, {l (/i^^) : < e < R} has the distribution of a Brownian 
motion up to a non-random time change. 

Second, the formulas (|2.9p and (|2.10p indicate that the covariance function does 
not depend on £2 when B^^^ (y) C B^^ (x) and does not even depend on £1 or £2 when 
(y) n (x) = 0. If one does the calculation with circular averages of the two 
dimensional GFF in each case analogously, one would see the same phenomenon. 
Namely, the smaller radius does not appear in the covariance function if one circle is 
entirely contained in the disk bounded by the other circle, while neither of the radii 
matters if the two disks bounded by the two circles respectively don't intersect. 
Such properties are consequences of the mean value theorem applied to the Green's 
function of the Laplace operator A in two dimensions, and these special properties 
of harmonic functions are no longer available to us in four dimensions. Nonetheless, 
we have seen from the above that by substituting {l (/i^ij) : x G M^, £ > O} for the 
circular average, we will recover in the four-dimensional setting properties similar 
to those in two dimensions. 

Finally, by examining the asymptotics of the Bessel functions near the origin, 
one finds that function G as defined in p.7p is asymptotic to — 2^ log r when r is 
small. Therefore, in both case (1) and case (2) from above, the covariance function 
is asymptotic to — log £1 for sufficiently small £1, while in case (3), the right hand 
side of p.lOp is asymptotic to — ^'^S \ ^ ^ vl for sufficiently small |a; — yj. In this 
sense, the covariance function of {l(/i^x) : x g R^,£ > O} does have logarithmic 
growth near diagonal as one would have expected. 

By now, one should be able to believe that {l(/i^x) : x G R"',£ > 0} is a rea- 
sonable replacement for the circular average of the 2D GFF in order to construct 
a random measure on M^. Indeed, the construction based on this Gaussian family 
will be carried out in the next section. We close this section with an important 
observation about T(/i^x) as a mapping from the variable x € R'^ to a random 
variable on Q under W. 

Corollary 3. Given £ > 0, the mapping a; € 1— )■ I {h^^) G (W) is continuous. 
Moreover, if a £ (O, |), then for almost every 9gQ,xE M.'^i-^I (/i^j) (9) £ R is 
Holder continuous with Holder constant a. 
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Proof. By Kolmogorov's continuity criterion ( ^Sllj . §4.3) applied to Gaussian ran- 
dom variables, in order to prove both statements in Corollary [3l it would be sufB- 
cient if we can show that there exist constant /? > and < C^.e < oo such that 
for every x,y G M^, 

To simplify the notations, we write as fif = fi (e) + /2 (e) da^, where both 
/i and /2 are actually analytic functions in e e [0, R] (the explicit formulas for /i 
and /2 are given by (17. 9p in the appendix). Therefore, we only need to show that 

both llcrf — and || dcr^ — c?crf||^_2 are bounded by Cp^^lx — yf ■ Perhaps 

the most straightforward way to see this is writing down the integral expressions 
for llcrf — erf ||^_2 and ||dcr^ — (io-f||^_2. Again, we refer to the formulas (|7.ip - (|7.6p 
in the appendix. From there, together with the series expression for the Bessel 
functions ([WatJ, §3.1), it is an easy matter to check that one can get the desired 
upper bound for \\a^ — (t^||^_2 and \\da^ — c?cr^||^_2 so long as /? e (0, 1). □ 

3. Construction of Random Measure 

In this section, we will use the Gaussian family {l(/i^j) : x e M'*,e > O} to 
construct a random measure on which "formally" can be written as "m (dx) = 
g27S(a;)gj2;" where 6* € 8 is chosen under the distribution of W, 7 > is a constant 
and dx is the Lebesgue measure on R*. Recall that at every x S /if tends 
to 2Sa: as e I in the sense of tempered distribution, so we can take the value of 
the random variable 9 ^X(ft,^j) {9) as the approximation for "9 (x)" as e J, 0. 
Furthermore, Corollary [3] certainly guarantees that with e fixed, one can always 
make the mapping 

{x,9) e R" X e^l{hf,.) (9) e R 

measurable with respect to *8k4 x *Be- In addition, we can assume for every 9 (z 
the function 

xeR^^ El [x) = exp {^X {h^.) (9) - ^—G (e)^ 

is positive and continuous, and hence if m\ (dx) = [x) dx, then {dx) is a 
positive regular and cr— finite Borel measure on R^. Moreover, given any Borel 
SCR"*, the mapping 

6* e e (B) [ El {x) da; e R 

J B 

is also non-negative and measurable. Hence by Tonelli's Theorem and the fact that 
[El {x)] = 1 for every x € R" and e > 0, 

[ml [B)] = /" [El (x)] dx = vol {B) . 

J B 

Since ml {dx) is simply the Lebesgue measure on R'' when 7 = 0, from now on 
we will only consider the case when 7 > 0. We want to study the convergence of 
ml {dx) as e I 0. In order to have the desired convergence, we only consider e taking 
values along a sequence {e„ = Cq : n > 1} for some fixed eo G (0, 1). Without loss 
of generality, we will assume ml^ {dx) is well defined as above for all n > 1 and 
every 9 £ Q. For the sake of convenience in later discussions, we will abuse the 
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notation by identifying "TOcq (da;)" with 0. We want to show that as n — > oo, al- 
most surely the sequence {m^^ {dx)} converges weakly to a limit measure {dx), 
written as m^^ (dx) mP (dx), in the sense that J^^ f (x) m^^ (dx) converges to 
/r4 / (x) 'mP {dx) for every f G Cc (M."^) which is the space of continuous and com- 
pactly supported functions on Mf^. To reach this goal, it suffices to show the con- 
vergence of /jj4 / (x) m^^ (dx) when / is any continuous function on F for any given 
compact set F C R''. In fact, we have the following result that holds for more 
general / so long as / is bounded and measurable on F. 

Lemma 4. Assume < 7^ < 2tt'^ and F C is compact. There exists a square 
integrable random variable G 1— > mP (F) G M+ such that 

N-l 

jm^^^j^ (F) — m^^ (F) converges to (F) as N ^ 00 

almost surely as well as in (W). 

Similarly, for every bounded measurable function f with supp (/) C F, there 
exists Af* (/) G (W) such that 



(/) = f (x) m^^ (dx) converges to M (/) as n —> cxi 
almost surely as well as in L? (W), and (/)| < ra^ (F) ||/||„ almost surely. 

Proof. To prove the first statement, we rewrite (^m.e„+i (r) ~ 'm-l^ (r)j , n > 1, as 
the following double integral: 



El^^ (x) El^^ (y) + El ix) < (y) - 2<^^ (x) < (y) dxdy. 

(pa L -I 

By Tonelli's Theorem, the W-expectation of above equals 
We split this integral by dividing the domain into two parts: 

// ''^'"^ //o<|x-y|<2.„ • 

The formula (|2.10l) implies the integrand is always zero in the designated domain of 
the first part. As for the second part, the integrand is always bounded by 2e'*'^'-'^'"+i' 
while the volume of the integral domain is bounded by Ce^ for some constant^ C. 
Together with the observations made in Section 2 about the asymptotics of G, one 
finds that 



(3.1) 



ml^^ (F) - ml (F) 



< 



(7g-(8^'-7')G(e„)^ 



Throughout this section, "C" denotes a positive finite constant that is universal in e„. The 
value of C may change from line to line. 
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The square root of the right hand side of p.l|) is summable in n > 1 and meanwhile 
m^^ (r) is clearly square integrable, so 

oo 

(3.2) ^(r)^^|m^„^^(r)-m^jr) 



n=0 



is square integrable and the convergence takes place in (W)- Furthermore, the 

series X)^=o (^) ~ (^') converges to (F) almost surely along a 

subsequence, but since the series is monotonic in N, it must converge to (F) 
almost surely along the full sequence. 

The second statement of the lemma follows from the same arguments. In fact, 
given a bounded measurable function / with supp (/) C F, if one replaces m^^ (F) by 

M^^ (/) in every step of the proof above, one can see that X^^^Lo -^e^„+i (/) ~ -^^cL (/) 
is also square integrable. The rest of the proof is straightforward. □ 

Theorem 5. Assume < 7^ < 2tt^. For almost every 9 E Q, there exists a 
non-negative regular a— finite Borel measure (dx) on R"' such that 

m^^ (dx) (dx) as n 00, 

and for every compact set F C , H^Ti^H^^j^p < (F) where ||'Ti^||^^j,p is the 

total variation of {dx) over F and nrfi (F) is as defined in ^3.2]) . 
In particular, for every f £ Cc (K^), 

/ (x) (dx) converges to f (x) (dx) as n 00 

almost surely as well as in (W). 

Proof. Clearly we only need to prove the first statement of the theorem, because 
assuming m^^ (dx) (dx) almost surely, the second statement is simply repeat- 

ing the second result in Lemma g] with (/) = J^^ f (x) (dx) for / G Cc (M^). 
As we mentioned earlier, to obtain the limit measure (dx), it suffices to show 
the convergence of mf^ (dx) on any compact set F C M^. We will achieve this goal 
via the Riesz representation theorem. We have already seen from the second part 
of Lemma m that, if denote M^^ (/) = J^^ f (x) m^^ {dx) for every n> 1 and every 
bounded measurable function / supported on F, then 

(3.3) M" (/) = lim Aff (/) exists and \m'^ (f)\ < ||/|l ^(F) < 00 

with probability one. However, to get the almost sure existence of (dx), we need 
to argue that with probability one, the statement above holds simultaneously for all 
functions / from a "large enough" class. To this end, we make use of the separability 
of the Banach space C (F) (when equipped with the uniform norm ||-||„). Namely, 
we can choose a countable sequence {fk : fc > 1} which is a dense subset of C (F), 
so almost surely the statement p.3p hold^ simultaneously for all fk, k > I. 

Now let's focus on 6* € 8 such that (|3.3p holds for all fk, k > 1. Given a general 
/ e C (F), let {fkj : J > 1} be a subsequence such that fk^ f in \\-\\^ as j — ^ 00, 



■''In this discussion, we will simply assume / = outside F for every f € C (F), so (/) and 
Mf (/) are still well defined. 
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then for every l,n> 1, 

\Ml{f)-Ml{f)\ < |Mf,(/)-A/f,(/,^.)l + K(A-,)-<(A.)l 

+ |Mf„(/feJ-Mfj/)| 
< 2^ (r) 11/ - f,^ \l + \Ml - Ml (/fej I . 

Obviously {M^ (/) : n > 1} forms a Cauchy sequence which immediately implies 
that (|3.3p also holds for / and in addition (/) = linij^oo {fkj)- Therefore, 
for almost every 9 G Q, f i-^ (/) is a linear and bounded functional on C(r), 
which, by the Riesz representation theorem, gives rise to a unique regular Borel 
measure (dx) on T such that (/) = )p / (x) (dx) for all / e C (F) and 
the total variation of (dx) is equal to the operator norm of which is bounded 
by (r). It's also clear that (dx) is non-negative and m^^ (dx) {dx). □ 

Compared with the second statement in Lemma 21 the second statement in The- 
orem[S]seems "weaker" since we have restricted ourselves to f G Cc (Mf^). However, 
we point out that the same statement, i.e., (/) converges to / (x) mP (dx) 
both almost surely and in (W) , is no longer true if / is only assumed to be 
a bounded measurable function with compact support. The reason is the follow- 
ing: for such a function /, although the existence of (/) = lim„^oo AI^ (/) 
is guaranteed by Lemma |31 in general one cannot draw any conclusion on the re- 
lation between (/) and J^^^ f {x) mP {dx), because m^^ {dx) only converges to 
mP {dx) weakly and one does not have control over — supp(/)' 

ever, under some circumstances, we can derive a relation between the two random 
variables {f) and f {x) mP {dx) . For example, if / = xa is the indicator 
function of a bounded open set A C K'^, then the weak convergence result implies 
{A) < lim„_j.oo m-e,, {A) almost surely. Meanwhile, the (W) convergence of 
m^^ {A) certainly leads to 



E 



w 



hm {A) = lim E^ [m^ {A)] = vol {A) ; 



on the other hand, let 1} C Cc (R"') be a sequence such that < fi XA 

as / — >■ cx), then by the monotone convergence theorem and the second statement 
in Theorem [5l 



E^ [m" {A)] = lim E^ 

l—¥00 



fi {x) m" {dx) 



i- 



= lim / /; {x) dx = vol {A) 



This can only be possible if {A) = lini„^oo {^) almost surely. More gen- 
erally ( |S11| . §9.1), if / is bounded and upper semicontinuous (or lower semicon- 
tinuous or m*— almost surely continuous), then it follows from a similar argument 
that (/) — lim„^oo Aff (/) almost surely, so the second statement in Theorem 
[S]also holds for /. 

The fact as stated above that E^ [m^ {A)] = vol {A) for every bounded open 
set guarantees that the limit measure {dx) cannot be almost everywhere triv- 
ial, i.e., W [m^ {dx) = O) < 1. In fact, we will prove later (in Lemma [TU)l that 
W {mP {dx) = O) =0, so mP {dx) is almost surely a positive measure. On the 
other hand, the following simple observation shows that mP {dx) almost surely does 
not assign positive mass to any given point. To see this, recall the assumption 
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< 7^ < 27r^ and the fact that G [r) is asymptotic to — log r when r is small. 
Then it is an easy matter to check that for any fixed a; G M"', 



E 



w 



0. 



limsupe^T'^^^-^m^ {x)) 
Therefore, if denote 

(3.4) = e e : Jim^ e^^'^^^-'m^ [b,^ (x)) = o} , 

then Qx is clearly a measurable subset of 8 and W {Qx) — 1- 

We will close this section by a remark about the condition of the constanl]^ 7. 
Readers may have noticed that the constraint < 7^ < 27r^ in Lemma 2] and 
Theorem [S] is more than what the proofs require. However, one reason of having 
this condition on 7 is that it guarantees the proof of Lemma |4] being correct even 
if one replaces 7 by 27. In other words, if we denote 



and define m^-^'*' (F) similarly using m^'^'^ (T) for any compact set F C R'', then 
77iS,27 (-p-j g^jgQ square integrable and in particular to^-^'*' (F) is almost surely finite. 
Some proofs in Section 5 make use of this consideration and hence the condition 
< 7^ < 27r^ becomes necessary there. We will remind readers when it comes to 
those situations. 

4. KPZ Relation 

Throughout later discussions, we will always assume < 7^ < 27r^ and for every 
9 d Q, (dx) is a non-negative regular and cr— finite Borel measure on R"* and 
m^^ (dx) (dx) (otherwise one simply assigns m^^ (dx) = mP (dx) = dx for all 

n > 1 on a measurable null set of 0). In this section, we would like to establish 
a KPZ relation between the Euclidean scaling exponent of a bounded (fractal) 
Borel set on K'* and its counterpart under the random measure (dx). We first 
recall from |DS11| some definitions. Given a bounded Borel _D C R"', the constant 
K E [0, 1] is called the Euclidean scaling exponent of D if 

log vol (Ja) _ 

illii . — rh. 

AiO log A* 

where A > and Dx = \Jx^dB\ (x) is the canonical A—neighborhood of D. In 
the random measure setting, for every A > 0, we first consider the mapping from 

r4 X e to [0,00]: 

\ sup {r > : [Br (x)) < A} if 6* e 6^, 
1 otherwise. 



(4.1) {x,0)^rA{x,e) = 



where Qx is as determined in (|3.4p . then define the isothermal A— neighborhood of 
i^by 

f either r\ {x, 9) > and dist (x, D) < r\ {x, 0) ~| 

(4.2) Z?'=<a;GM: >. 
^ ' I or TA (x, 6*) = and X e D | 



"In this article, we don't have particular emphasis on the potential physics meaning of 7. 
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Intuitively, B^f^^xfi) (x) is the "largest" ball (in radius) centered at x with volume A 
under the random measure (dx), and D^'^ is the "neighborhood" obtained by 
covering D with all such balls with equal volume. If there exists constant K e [0, 1] 
such that 

logE^ [m« (D^'<')] 

4.3 lim ^— ^ ^ = K, 

AiO log A 

then K is called the quantum scaling exponent of D. K can be viewed as the 
"expected" scaling exponent and as such the counterpart of k under the random 
measure mP (dx). In two dimensions, k and K satisfy the so-called KPZ relation 
which is a quadratic relation. Our goal in this section is to extend such a relation 
to the four dimensional setting. 

However, we haven't yet justified the "meaning" of the notations in (|4.3p . First 
we have to check that D^'^ is a Borel set in R"* for every 6 £ Q, which requires 
us to verify the measurability of {x,0) H' rA{x,9) with respect to ^Br4 x Se- 
To this end, we observe that {x,9) i— )■ {B^ (x)) is measurable for every r > 
because there certainly exists continuous mapping x £ i-^ pf £ Cc (W^) for 
every I > 1 with < pf Xb,.(x) and hence (pf ) /• {Br (x)) as ^ — >• oo 
for every {x,9). From this we conclude that {{x,9) : 9 £ Qx} is a Borel set in 
M'* X 0, and then the measurability of ta {x, 9) follows simply by identifying the set 
{{x, 0) : < TA (x, 9) < a} with {{x,9) -.9 £ <dx,m^ [Ba (x)) > A} for every a > 0. 
Therefore, for every 9, D^'^ is a Borel set in since it is the 6*— section of the 
following Borel set in R'' x 0: 

(£» X 9) U {{x, 9) : rA {x, 9) > 0, dist {x, D) < ta {x, 9)} . 

Next, we need to show that 9 i-> (Z)'^'^) is measurable with respect to Se- 
More generally, if C G *8r4 x So and ee G R"* : (x, 9) £ C} is the ^-section 
of C, we claim that 9 i-> {^^) measurable with respect to !Be- To see this, 
denote = H Bn (0) for every iV > 1 and choose a sequence of mollifiers 
{(pk £ [0, 1] : fc > 1} C Cc (W^) such that = tpk * Xc% converges to Xc% under 
\\-\\^ as fc — > oo. Moreover, for every fc > 1, since C £ *Br4 x 58e, 

J {y.e)ec,\y\<N 

is also measurable with respect to *Br4 x *Be- Thus, because g^ £ Cc (W^) and 
< 1 for every fc > 1 and m^^ (dx) mP [dx), 

(C") = lim lim hm / gl (x) E^^ (x) dx, 

N—>oo fc— >oo n— >oo Jjj4 " 

which is a measurable function in 9. 

At this point, one can already tell that it is convenient to consider the spa- 
tial variable x and the GFF 9 at the same time. In particular, if Ai {dxd9) = 
(dx) W {d9), then it's clear from the preceding that A4 {dxd9) is a non- negative 
(T— finite Borel measure on the product space R'' x Q. To connect k with K, it 
is natural to investigate the distribution of ta {x,9) under the M. {dxd9). To get 
started, we will first look at the distribution of {Br {x)) under {dxd9) for any 
given r > 0. The following lemma says that we can change our "perspective" by first 
choosing a "base" point x £ Mf^ and then examine the distribution of mP {Br {x)) 
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at X. This is reaUzed by a procedure of exchanging the order of integration. The 
proof of Lemma [6] is given in Section 5. 

Lemma 6. For every x € M^, let be the measurable subset of Q as determined 
in ^3-4^ , and define 

'exp(^^Ko{\x-y\))m'idy) tfOeQ,, 
^m^ (dy) otherwise. 



(4.4) e ^ m**'^ (dy) 



Then almost surely rh {dy) is also a non-negative regular and a— finite Borel 

p4 



measure on . 



Moreover, for every r > 0, compact set F C and F £ Cq (M!^ x [0,oo)), 
(4.5) 

[ [ F {x,m^ {Br{x)))m'^ {dx)W {d0) ^ [ [ F {x.rh^''' {Br {x)))yV [dO) dx. 
Je Jr Jr Je 

In particular, this implies that for every r > 0, the joint distribution of {x, {Br {x))) 
under M {dxdO) is the same as the distribution of {x, rh • ^ {Br {x))) underW{d9) dx, 
whose marginal distribution on Q at x is independent of x. 

Based on Lemma [6l instead of {Br {x)) under {dxdO), we may as weU 
study the distribution of 7n^'^ {Br {x)) under W {dO) dx. Similarly, to understand 
ta {x, 9) under M {dxdO), we only need to look at the random variable given by 

, m ~ / m fsup{r >0:m'''-(B,(x)) < A} if G 6,, 

(4.6) {x,9)^rA{x,9)= \ ^ ■ 

I otherwise. 

under W {d9) dx, whose marginal distribution on O at a; is again independent of x. 

To proceed from here, we will follow the same strategy as in [DSllj . For the 
sake of completeness, we will still present the main steps here. For every r > 
and A > 0, since the distribution of rh^'^ {Br {x)) and ta {x,9) under W does not 
depend on x, we can assume x is the origin and simplify the notation by denoting 
Br = Br (0), f-A {9) = fA {0,9) and {Br) = rh'^^° {Br (0)). We want to find an 
"approximation" for rh^ {Br) by conditioning on the value of the GFF restricted 
to the "boundary" of Br. To be precise, recall from the definition (|2.1ip that if 

r {t) = {t + G {R)) with t > 0, then Xt = I (^/."^^J - ^ ^^s the same 

distribution of a standard Brownian motion. We want to investigate the conditional 
expectation of {Br(t)) given Xt. To do this, we need to relate mP {^Br(t)) to the 
approximating measures to^^ {dx) , which requires us to overcome the singularity of 

62^^°"'''' at the origin. To this end, assume {/;:/> 1} C Cc {Br) is a sequence 
with < fi X-Br(t) '"^^ ^ therQ 

a (S^) 3 d, {■) ^ f, {■) ^ ^^^^^^ {.) ei^^od-l). 

Therefore, one can apply the convergence results in Theorem[S]to di for every / > 1. 
Together with the monotone convergence theorem, one sees that for every t > 0, 

(4.7) [m^ {Brit)) \Xt] = lim liminf [Mf (d,) \Xt\ . 



■^The notation "a V /3" denotes "max {a, /3}" and "a A /3" denotes "min {a, /3}". 
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Given the order of taking limits in the right hand side of ()4.7|) . for every I > 1 
and eventually all large n, we have B^^ (y) C 5^(4-) C Bn for every y € supp (d/). 
Thus by a simple exercise on conditional expectations of Gaussian random variables 
along with (|2.9p . we can derive from (14.71) that 

e^Ko(\y\) [(^^ (Ij^i) _ p 



(4.8) 



• exp 



^^(/o(|y|)-/2(|y|)PWf 



dy, 



where P it) = (AttH) ' [{if - loh) \r {t)) - {if ^ hh) ' (R) . 

If one carefully examines the asymptotics near the origin of the Bessel functions 
involved, one realizes that (|4.8p suggests the conditional expectation of rh^ 
given Xf, when t is large, is "approximately" 

(4.9) ■ ^ / - -2 



m'^* (Brit)) = exp jX, 



For the moment we will "pretend" rh^ (^Br(^t)} is just m 
KPZ relation under this circumstance. 

For every A > 0, we define the stopping time: 



e* 



Br{t)) and formulate the 



(4.10) 



Tt 



inf U > : m®* {B, 



r{t)} 



exp ^X; 



St:' 



T 
2 



t] <A 



The distribution of can be completely determined by a martingale argument, 
which is straightforward but worth repeating. Namely, for every s < 0, by Doob's 
stopping time theorem, |exp 



sXfATi — ^ (t A T^) : t > o| is a uniformly bounded 
martingale. Furthermore, the continuity of Brownian motion implies that 

log A o Tt 



1 



2 



_A 
7 ' 



Therefore, the fact that the expectation of the martingale at t = is equal to that 
at t = leads to the formula of the moment generating function of T^: 



(4.11) 



E 



w 



7s 



exp -- 



2s Svr^ - 



27 



Tt 



From here we can derive our first version of the KPZ relation which is easy but 
revealing. 



Lemma 7. Assume D C 
K e [0, 1], i.e., 

(4.12) 



is a bounded Borel set with Euclidean scaling exponent 

limi^^i^... 

a;o log A* 



For every A > 0, let T^ be as in l{4.10^ and define the random "radius": 

e^rl {e)^G-'{Tl{6) + G{K)), 
and the random "neighborhood": 
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The 



have 



lim^ ^=lini , . 

A^O log A AiO log A 



(-a)^ 



(4.13) lim^ k e ^=lini ^-^ = K 

where K G [0, 1] is determined by the following quadratic relation with k: 



(4.14) 



K 1 



7 



167r^ 



7 



167r2 



Proof. Clearly in this setting, we want to cover D with open balls that all have the 
same "critical" radius determined by the stopping time associated with the martin- 
gale rh^* (5^(4)) (defined in (|4.9p '). We don't have to worry about the measurability 
of 61 (D^*-^) because both ^ rl [6) and (6',r) ^ mf> (Dr) are measurable. 

Conditioning on r]!^ = r G (0, i?], D^*'^ is bounded and open, and the conditional 
expectation of (^D^*'^) is a multiple (reciprocal of the probability density func- 
tion of at r) of vol (-Dr), which, according to (|4.12p .R'?^" r'^" for every r G (0, R]. 
This means E^^ [m^ (D^*-^)] « E^^ Urlf"] and further w E^ [exp (-Stt^kT^)]. 

Given (|4.1ip . clearly one wants to set Stt^k to be ^ ^ (stt^ ^ 'hr) some 

s G 



"7, 0], in which case 
logE>^ [m<' 



lim ■ 

AiO 



log A aJS log A 

The results in (|4.13p and (|4.14p follow immediately after setting K = —-. 



logE^ [exp (-Stt^kT;)] 



s 



□ 



Next, we argue that -mf* {Br) is indeed a "legitimate" approximation for rh^ {Br) 
in the sense that r\ , as defined in Lemma [71 approximates rA-d'^rA {S) when 
"compared" in the limit of the logarithm ratio. 



Lemma 8. Assume the pair (k, K) 
(TW . Then, 



G [0, 1]^ satisfies the quadratic relation in 



(4.15) 



lim ■ 

A4.0 



logE^ ( 



rA) 



log A 



K or equivalently lim ■ 



log 


;EW 


'{fA)'\ 


log 


;EW 


\rl)'\ 



The proof of this lemma is given in Section 5. There we also prove a preliminary 
result (Lemma [TU)) which actually implies the almost sure non-triviality of the 
measure mP {dx) as well as {dx) . Most importantly, this lemma builds up the 
final passage to the KPZ relation for {dx), the "true" case in which we are 
interested. Again, we will only present the statement here and leave the proof to 
the next section. 

Theorem 9. Let Z? C R** 6e a hounded Borel set with Euclidean scaling exponent 
K G [0,1]. Then D has quantum scaling exponent K G [0,1] as defined in \4-^ ! 
where K is related to n by ^.14\ l- 



Throughout this article, the notation "ft:" means "bounded from above and below by a uni- 
versal constant multiple of". 
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5. Proofs of Results in Section 4 

We will now prove Lemma [6l The strategy is to relate {dx) to the approx- 
imating measures m^^ (dx) and recognize that, by the Cameron-Martin formula, 

the density of m^^ (dx), i.e., E^^ (x) — exp (^jX (^^/^j ^ (^) ^ (e„)^ , is just the 
Radon-Nikodym derivative with respect to VV of the Gaussian measure induced by 
the translation 6 9 + under W. Note that the constraint < 7^ < 27r^ 

becomes necessary in this proof. 

Proof of Lemma\^ Let mP'^ (dy) be the measure as defined in (j4.4l) . The claim 
that rh^'^ (dy) is non-negative regular and ct— finite follows from the observation 

that exp {^-^Kq {\x — -j)^ is locally integrable with respect to mP (dx) if e 6^;. 
Without loss of generality, we will assume x = 0. The only possible problem comes 
from the singularity at 0. However, if we rewrite 

2 00 „ 2 

B.o(O) fc=o-^«fc<lal<e(c-l 

00 2 

<Y^,^Koi^^)m<^ (i?e,_,(0)), 
fc=0 

then the criterion p.4[) for Qx guarantees that the series in the right hand side of 
above is convergent. 

Now we move on to the second part of the lemma. Clearly both mappings 

(x, 9)^ F (x, {Br {x))) and (x, 9) ^ F (x, 771**'"= {Br {x))) 

are measurable with respect to *Br4 x Se, so the two integrals in (|4.5p are well 
defined and in fact finite. Choose a continuous mapping x £ F G Co (M*) with 

< < XBr(2;)- We first show (|4.5p holds with x_b^(2;) replaced by p^. Namely, 
we claim that 



(5.1) 



F (x, (p^)) (dx) W {d9) = 

F (^x,Af« (^p^ei^^^fl^^-l)^^ W(d6i)dx. 



r Je 



We start with rewriting the left hand side of (|5.ip . Since F{x,AI^{p^)) is 
continuous in x G F, the weak convergence result implies that 



F (x, (p^)) (dx) = Mm F (x, Af® (p^)) {dx) , 
which, by the dominated convergence theorem, leads to 
F (x, (p^)) (dx) W (d^) 

(5.2) 



e 



n— J-oo 



lim / / F [x^Kf {p'))El {x)dxW{d9) 



By Fubini's Theorem and the consideration (about viewing (x) as the Radon- 
Nikodym derivative of the translated Wiener measure) we made before the proof. 
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we have that the right hand side of (|5.2p equals 



hm 



r 



(x,M 



Now given x G M'*, the Cameron-Martin theorem guarantees that also with proba- 
bility 1, rric^ (dy) weakly converges to m^^''''^«" (dy) as fc — > oo simultaneously 
for all n > 1. In particular, 



k^oc 



= lim 



(y) exp b^E 



I h 



I h 



El{y)dy. 



At this point, it is clear that (jS.ip would follow if we can show that for every 9 € O^;, 



(5.3) 



lim lim / p'^ [y) exp {-f^'R \x (h^y )l(h 

n— >oo fe— >oo \ I \ '^'k J \ ' 



El{y)dy 



j^^ p"" (y) exp {^^Kn {\x - {dy) < oo. 



The right hand side of (|5.3p is finite for e as we have seen in the proof of 
the first part of this lemma. To establish the equation in (j5.3p . we assume n > 1 
is sufficiently large and k > n and divide the integral in the left hand side of ([57 
into three pieces: 

(5.4) 



j/-x|<e„-efc J €„-ek<\y-x\<en+ek J \y-x\>e„+ek ) 

p- (y) exp Ll (^h^,y^ ) (9) + ^^E [l (^h^y^ ) I (^h^.^ ) 



YG(efc)|dy. 

We will investigate the limit as fc — >■ oo and then n — >■ oo of each piece separately. 
By (|2.10p . the last integral in equals 

(5.5) / p-(2/)e^^°(l^-^l)<(y)d2/. 

If the domain in (|5.5p is replaced by {y : |j/ — x| > e„}, then the integral would be 

{y)e^''^'^\--y\^El{v)dy = 



(5.6) 



\v-x\>e 



p- {y)El{y)dy. 

As fc — > oo and then ?i — > oo, the first term in the right hand side of (|5.6p converges 

2 

to p^ {y) ^°'^^^~y^'' mP {dy) (which is the term we want and also the only 
term that should survive in the limit). On the other hand, as fc — > oo, the second 

term on the right hand side of (|5.6p is bounded by e^^°^'^'^'^mP [b^^ {x)^ , which, 
because 9 £ converges to zero when n — >■ oo. As for the "redundant annulus" 
which is the difference between the left hand side of (|5.6p and (|5.5p . it's bounded 
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by e2^^°*-^"-* times the integral of E^^ over the annulus {e„ < |a; — -j < e„ + 6^}. 
One can apply the Schwarz inequality to see that this integral is bounded by 

(5-7) e^^(^^)vol({e„ < \x - ■{ < en + ek})^ (^^(2r))' 

where 2r = {2y : y G F}. Given the considerations at the end of Section 2, without 
loss of generality, we can assume m^'^i (2r) is finite and hence the limit of (|5.7p as 
fc — > cx) with n fixed is zero (because the volume of the annulus « « e^^'^ G(ek) 
and 27r^ > 7^), so this "annulus" is negligible. 

The second integral in (|5.4p becomes negligible by a similar argument. This time 
the Schwarz inequality implies that the second integral is bounded by 

g7VG(^'=)G(^")e^G(^.)vol5 ({e„ _ < - .| < e„ + e J) (^^^ (2r)) ' 

where we use the simple estimate E I {h^v^I {h^^ ^ < ^JlT(ekyG{€n) . Again, 
with n fixed, the factor that involves k converges to zero. 

As for the first integral, because of (|2.9I) and the asymptotics of the Bessel 

functions involved, E I {h^y ^ I (jif^^ ^ is bounded by r/G (e„) for some constant 
rj e (1, 2) for all sufficiently large n. Therefore, with n fixed, the integral as /c cx) 
is bounded by e^'^'^^'^^'^m^ (^e„ i^)^ j ^^nd as n — > 00 it therefore converges to zero. 

So far we have proved the claim (|5.1I) . To reach (|4.5p . one takes a sequence 
{pf : 1>1} C (R4) such that < pf XbAx) as ^ -> 00, and for every I > 1, 
X G i-^- pf e Co (R'') is continuous. So ()5.ip holds for each / > 1. After carefully 
examining the integrals on both sides of (|5.ip . one realizes that the limit as Z — ^ cxd 
can be passed all the way inside to produce (|4.5p . At the end, it's clear that given x, 
the distribution of m^'^ {Br (x)) under W is independent of x due to the translation 
invariance of measure (dy). Hence we have completed the proof to Lemma[B] □ 

Now we move on to the proofs of the KPZ results. The techniques we adopt 
here differ from those used in the two dimensional proofs in |DS11| , partly because 
of the absence in our setting of the two dimensional conformal structure as well 
as the compactness of the domain. For example, the next lemma is the "tail esti- 
mate" which is the key estimate in proving both Lemma [8] and Theorem [9l In the 
two dimensional counterpart, the corresponding estimate ([DSll], §4.3) is a super- 
exponential type of estimate. Below we prove an exponential type of estimate in 
the four dimensional setting, but by carefully "tuning" the exponential decay rate, 
we can still make it sufficient for our purposes. Again, the occurrence of (dy) 

in the next lemma refers to the measure e2^^''*-'^'^m^ (dy) assuming 9 & &o- In 
other words, only balls centered at the origin are concerned. However, since the 
distribution of ?7i^'^ (Br (x)) under W does not depend on x, the same result will 

hold for m^'^ (dy) = e^^^""^"^"™^ (dy) (assuming 9 e Qx) no matter what x is. 

Lemma 10. Let B be the closed ball in centered at the origin with unit volume 

2 2 

under e^^°^^^'^^ dy ,i.e., Jg e^^^^^'^'-'dy = 1. If 5 and p are constants satisfying 
< (5 < 47r2 - 272 and — ^ <P<1, 
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then there exists C > such tha^ for all sufficiently large A > 0, 



6 A 



(5.8) W (m" (B) < e-^^) <Cexp\-^( Stt^ _ ^2 _ 7 

L 7 V / 

Proof. Since B is closed, it suffices to estimate W (lini sup^^^^ {B) < e^^^) 
wliere rh^ (dy) has density es^^"'-'^''' with respect to (c'y)- I3y the same 



argument as used in deriving the estimate p.ip . we can show that there exists 
constant C > such that for all n > 1, 



E 



w 



< 



For any S with < 6 < Aii^ — 27^, denote A'^ n>l, the measurable set 



yi > n, 



(B) - ml (B) 



Then it follows easily from Chebyshev's inequality and the Borel-Cantelli Lemma 
that W (Ur=i -^n) = 1- Moreover, if Ai = A[ and An = A'„\A'„_i for n>2, then 
there exists constant C > such that for all n > 2, 

(5.9) W {An) < Ce2^7e-(8-^-7^-5)G(.„)_ 

Set B = {limsup„^^m«^ (B) < e-^T}, then W (S) = W ^ y^„) and 

it's clear that 9 g BnAn implies rh^ {B) < cse~^^ where cs = l + e~5G(cii)^ 
Given any p such that g^^_^^_^ < p < 1 (notice that such p always exists since 
< (5 < 47r2 - 27^ and 47r2 + 7^ < Stt^ - 7^ - 5), we set up the "threshold" N eN 
which is the unique (recall that G is strictly decreasing on (0, 00)) integer such that 



(5.10) 



G(e^)<^butG(eA,+i)>^. 

7 7 



The desired estimate (|5.8p is trivial when n > + 1, because (|5.9p and (|5.10p 
implies 



J2 WiBnAn) < Ce2'47e-(8-^-7=-5)G(e„ + J 



n=N+l 



< Cexp 



7 



Jtt^ - 72 - 2_ _ 5 M 



When n — 1, ■ ■ ■ , N, we apply Jensen's inequality to see that 
W (m^^ (S) < c^e-^^) 



(5.11) ^ ^ (,^^P 

< 



7X(/i^.J(0)--^G(e„))e- 



ri^o(lal) 



W (/^I (V„) (^)ei^-°(l^l)ci2/ < -A+lGien) + 

By Corollary [31 without loss of generality, we can assume that for all n > 1 and 
every 6, the function y B ^ I {h^i ^ [d] is continuous and hence uniformly 



Throughout this section, C denotes a constant that may depend on 7, (5, p and R, but universal 
in A, tn, X and A. The values of C may change from line to line. 
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continuous on B. Therefore, one can easily check (for example, by writing the 
integral as the limit of a discrete sum of Gaussian random variables) that 



dye 



is also a centered Gaussian random variable for every n > 1, and furthermore, the 
variance can be bounded by a constant M that is universal in n > 1. In fact, M 
can be taken as a constant multiple of 



Ko {\x - y\) e^(^"(l"l)+^«(l^l»da;dy. 



BxB 



Since G (e„) < ^ for n = 1, • • • , iV, A - (e„) > (1 - p)A, and (ICTD implies 
that when A is sufficiently large, 

W {ml^ (B) < cse"^^) < exp 



1 

<exp -^[{l- p)A- ^ Xogcs 



<exp pf A^ 



■ log cs 



In addition, (|5.10p implies that N is approximately a constant multiple of A. There- 
fore, when A is large, 



AT 



N 



^ W (S n X) < ^ W (m^^ (B) < cse-^^) < CAe 



n=l 



n=l 



So W (S n An) actually decays super-exponentially fast as A — )■ oo, and this 

estimate can certainly be transformed into the desired form as in (15. 8p . □ 



We are now ready to prove Lemma [S] Recall the notation ta : 6' (0, 0) 

where f\ (0, 6) is as defined in (|4.6I) with x being the origin. Let {k,K) e [0,1] be 
a pair as in (|4.14p . In order to get (|4.15p , it suffices to show that 



(5.12) 



(rA) 



4k 



< c 



for some constant C > universal in A as A 4. 0. We will prove the existence of the 
upper bound and the lower bound in (|5.12p separately. 

Proof of the upper bound in 15.12\) : For notational convenience, we introduce 
the "stopping time" corresponding to f\, i.e., T\ = G {f\) — G{R). We want 
to show A-^E^ \ exp (— Stt^kTa)] is bounded from above uniformly in small A. 
It's clear, from (|5.8p and the fact that W (Qq) = 1, where 9o is as in p.4p . that 
Ta G (— G(i?) ,oo) almost surely. Let constant S and p be as in the statement of 
Lemma [TUl Set 



S = 



-log A 



2p (Sn 



87r2 - 72 



2p (, 



87r2 _ 72 _ :^ _ ^ 
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Then the expectation of exp (— Stt^kTa) can be written as 

[exp (-Stt^a^Ta) X{-GiR)<Tj.<S}] + [exp (-Stt^a^Ta) X{S<Ta<oo}] • 

In the first term, Ta < S* implies that the volume of the closed ball centered at 
the origin with radius r (S) = G^^ {S + G (R)) ( « exp (—271^5*)) is no greater 

than A under the measure rhP [dy)^ while this ball has volume sa [r [S)) 

(~ exp (— (Stt^ — 7^) S") ) under the measure e2^l^^''"*''-'(iy. However, by Lemma 
[TOl the probability of this event is bounded by 

Cexp (S.^ -,^~:^-s) - '-^^S 



7 \ P / V 7 7 

which, given this particular choice of S, is equal to a constant multiple of A^. 
Therefore, the first piece of integral causes no trouble. 

The second integral is bounded by e~^'^ '^^ . Hence we only need to check that 
A~^e~^'^ '^^ , or equivalently, exp (— iiTlog A — Stt^kS) stays bounded as A J, 0. In 
fact, we will show that for all possible values of (k, K) and all sufficiently small 
A > 0, Klog A + Stt^kS > 0, that is (assuming log A < 0), 



8^' - 7' 2p (8^2 _ ^2 _ £ _ 

To simplify the notations, let's write C, = 2p ^Stt^ — 7^ — — 5^- Recall from the 

statement of Lemma [TUl that < ^ < An"^ - and < p < 1, so ( > Stt^. 

If we express k in terms of K according to (|4.14p . then the statement in (|5.13p is 
equivalent to 

F (K) = -f^K^ + {16tt^ -1^-C)K-C<0 

for all possible values oi K G [0, 1]. However, this is clearly true since F is quadratic 
and (0) = -C < as weU as F (1) = levr^ - 2C < 0. □ 

Proof of the lower bound in i5.12\} : Recall that is the stopping time (as de- 
fined in (j4.10p ) associated with the "approximating" measure rh^*, and 

[exp (~8^2^r;)] = A^. 

We observe that E^ [exp (— Stt^kTa)] is greater than the integral of exp (— Stt^kTa) 
over the subset {Ta < Ta}, where the integrand is greater or equal to exp (— Stt^kT^) . 
Therefore, we have 

E^ [exp (-Stt^kTa)] > E^ [exp (Stt^kT^)] - exp (Stt^kT^) X{ta>tx} ■ 

It is clear that in order to get the desired lower bound, we need to find constant 
< c < 1 such that 



(5.14) A 



exp (-Stt^kTa) x{ta>ti} 



< c 



uniformly in small A. Conditioning on = T, Tj\ > T implies {Bj.{t)) > A 
and hence ^^;) > 1. By Chebyshev's inequality, other than a factor given 
by the probability density function of T^, the conditional probability of {Ta > T} 
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is bounded by the expectation of -rjz^^'^^^^ which, given the expression in (|4.8p 

(which is the conditional expectation of the numerator given the denominator), can 
be bounded by constant c € (0, 1) which is universal in A and T. So the estimate 
in (|5.14p will be satisfied by this choice of c. □ 



Proof of Theorem\^ Assume D C (0) for some sufhcently large > 1. Let 
TA (a;, e) and D'^'^ be as defined in gl]) and (|42]). Denote TV (^6*^2;) ee W {dO) dx. 
Based on LemmaH [m^ [D^^")] equals 

M {{{x, 9) : either x € D or dist {x, D) < ta {x, 9)}) 

(5.15) = Mi{ix,9):\x\<2N,dist{x,D) <fAix,9)}) 

+ lim M ({(x, 9):2N < \x\ < M, dist {x, D) < ta (x, 9)}) . 

In the first term of the right hand side of (|5.15p . conditioning on fjy {x,9) under 
M {d9dx), since its marginal distribution on Q does not depend on x, the conditional 

probability of the set is proportional to vol ^-DfA(2;.6i) H B2N (0)^ . We further split 
the set into two cases: ta (x, 9) > N and ta (x, 9) < N, the later of which also 
implies Dfj^(x,e) ^ B2N (0). Therefore, the first term can be rewritten as (up to a 
constant depending on A^) 

E'^ [vol (Df^i^x^g)) X{fA{x,e}<N}] + E^ vol (^Df^^^^g^ n B2N (0)) X{fAix,e}>N} 

[vol{Df^^a:,e))] -E'^ [vol(Df^(^,e))x{fA(a:,9)>Ar}] 



vol (^L'fA(2;,e) n B2N (0)^ X{fA{x,e}>N} 



According to the assumption (|4.12p and Lemma|51 E^ [vol (i^fAC^.e))] ~ when 
A is sufficiently small. On the other hand, given {x,9) > N, Df.A{x,9) is always 
contained in the ball centered at the origin with radius 2f\ {x,9), so the last two 
terms in the right hand side of the equation above are both bounded by (up to a 
constant) 



(5.16) E^ (fA(x,0))"x{fA(x,9)>A^} <4 / u^W{rAix,9)>u)du. 

L J ^[1,00) 

If C = 2/9 ^Stt^ ~ 7^ ~ — where S and p are the same as in the statement of 
Lemma [TUl then by (15. 



W(fA(x,6') >u) < W(m'''^(B„(a;)) < A) <CA^u 

Given the particular range of (5, p and one sees that not only is the integral in 
(|5.16p finite, but it also converges to zero faster than A^ as A J, for any possible 
value of a: e [0,1]. 

In the second term in (|5.15l) . since D C B^ (0), the assumptions \x\ > 2N and 
dist (a;, -D) < ta (a;, 0) imply fA{x,9) > ^\x\ whose probability, as we have seen 

earlier, is bounded by CA->^ \x\ t^V 2^ J which is integrable (with respect to dx) 
ill the entire domain {|x| > 2A^}. Therefore the second term also converges to zero 
faster than A^ as A | 0. To summarize, we have shown that [m^ (Z)^'^)] is 
a constant multiple of A-^ + o (A^) as A J, which is sufficient for the desired 
conclusion. □ 
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6. Possible Generalizations 

Generalizations to M^": In this subsection, we outline a possible generalization 
of the four dimensional treatments carried out in previous sections to higher even 
dimensions R^" with n > 2. We consider the GFF on R^" with the underlying 
Hilbert space H = i/" (R^") which is the completion of the Schwartz test function 
space S (R^") under the inner product ((/ — A)" •, •)^2 • Similarly, for every x S M^" 
and e > 0, af denotes the tempered distribution which is to take the spherical 
average of a test function over the sphere Se (x). In this setting, G H~" (M^") 
and again if ha-^ = (/ — A) "erf, then ha^ G H and the Paley- Wiener integral 
I [ha^) can be viewed as the "generalized" action of on the GFF. Moreover, 
the higher order of the operator (/ — A)" allows us to take higher "derivatives" of 
erf in the radial variable e. If (i"Vf = j-^cr^ is defined in the sense of tempered 
distribution for every m G N, then simple computations of the Fourier transforms 
show that 

(0 = C„e^(-'«)«- (6|^|)i-"j„_i(6|^|) 

alidad a J [^) ^^^aaU-^«e de™ 1^ (e J' 

where C„ > is a dimensional constant. In particular, we can write 

where f |Wat| . §3.31) (p (r) = '^^'ii''-' for r > and (^(™) (r) is analytic in r near 

and asymptotic to r as r — > oo for every to G N. Therefore, erf and d^crf 

for 1 < m < n - 1 are in i/"" (R^") . 

We can mimic the approach in Section 2 and define the vector- valued Gaussian 
random variable on 8: for every x G R^" and e > 0, 

It turns out that in this setting we can also compute the covariance matrix of 
the family {Vf : x G M^",e > O} explicitly under each circumstance as prescribed 
in Lemma [U and the covariance matrix also has a similar "separability" property 
as in four dimensional case. In fact, following a similar computation as the one 
(provided in the appendix) conducted to prove Lemma [TJ it is not hard to see 
that there exist invertible n x n matrices A (r) , B (r) , C (r) and D (r) for every 
r G (0,oo), such that all the entries of A (r) and D (r) are functions in the linear 
span of {r^'^Kk (r) : < i < k < 2n - 2} , while all the entries of B (r) and C (r) 
are in the linear span of {r-'^Ik (r) -.0 < £ < k <2n-2}. Moreover, given X (z K 

V,^ {VfJ^ = A(ei)B^ (ea); givenx,y G R^" with x 7^ y and 
[Ve, {y.lV] = A(ei)C(|a:-j/|)BT(e2); given x,y e M^" 
with X ^ y and \x-y\ > 61+63, [^^ (K^)^] = B (61) D (|x - y|) B^ (62). 

Therefore, if we similarly defined the "normalized" vector C/f = B^^ {e)V^, then 
the Gaussian family {/7f : a; G M'^",6 > O} will have the same properties as those 
of the corresponding family (also denoted by U^) in four dimensions. 

On the other hand, all the entries of the matrix B (e) are linear combinations of 
e^'/fc (e) with < I < k < 2n — 2, and if one lets 62 ^ in the covariance matrix 



and 61 > 62 > 0, E'^ 
61 > \x~y\+e2, E^ 
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obtained in the second circumstance (when ei > |a; — y\ +62) from above, combined 
with integral expressions for the entries of the covariance matrix, then one can easily 
conclude that there exists constant matrix B which is non-degenerate (hence so is 
B~^) such that B (e) converges to B as e 4. 0. Therefore, all the entries of B^^ (e) 
must be analytic in e near zero. In particular, by examining the asymptotics of the 
entries of B~^ (e) near zero, one can find the appropriate constant vector C e K^" 
such that (t/f , C)R2n "approximates" the GFF at x when e is small in the same sense 
as described in Section 2. 

Clearly, (C/f , C)r2" has all the properties of I (/i^^) in four dimensions as stated 
in Theorem H When ei > £2 > 0, 

when £1 > £2 + \x-yl [{Uf^X)^2„ (f^."^' Or^J = c(£i,|x-y|) which is in- 
dependent of £2; when |a; - > £1 + £2, C)r2„ (t^Jl' C)k2„] ^ d{\x- y\) 
which is independent of £1 and £2. In principle, we can derive the explicit formulas 
of G(£i), c{ei,\x — y\) and d{\x — y\), and one can expect that they have loga- 
rithmic growth when £1 and \x — y\ are small because the Green's function of the 
operator (/ — A)" on M^" has logarithmic growth near the diagonal. Therefore, 
it's reasonable to believe that if one takes (J7f,C)K2n to construct a sequence of 
approximating measures, i.e., 

ml {dx) ^ exp (^7 (f/f,, C)k.„ iO) - (£fe)) dx, 

then the sequence {ml : A: > l} will almost surely admit a limit measure in the 
sense of weak convergence. Furthermore, the quantum scaling component of a 
bounded set on M^" under this limit measure should also satisfy a quadratic relation 
with its counterpart under the Lebesgue measure. However, the amount and the 
complexity of computations quickly become considerable as n increases. 

Generalizations to Manifolds: In this last part we explain a more conceptual 
approach to constructing analogues of the two-dimensional GFF on compact even- 
dimensional manifolds. As we have remarked in the introduction, in dimension two, 
the GFF defines a measure on a conformal class of metrics on a Riemann surface E, 
constructed starting with a reference metric go on E, but in the end independent 
of go- In fact, the GFF inner product of two functions /i, /2 S C^(E) is defined 
by fpH lDSTTllDSOQllHMP] ) 

(/i,/2)a,„ ^ (/i, A90/2) ^ ljiix){Ag,h){x)dvo\g,{x), 

where Ag^ is the Laplace-Beltrami operator on E with respect to go. This inner 
product is conformally invariant. Indeed, if the metric go is changed conformally 
to gi = e^'^go for some uj £ (S), then the volume element changes as 

dvolgj = e^^dvolgg. 

while the Laplacian is changed as Ag^ — e^^^Ag,-, . Therefore, after obvious cancel- 
lations we find that 

It seems natural to define a similar measure for conformal classes of metrics 
in higher dimensions. Below, we explain how to do that for certain conformal 
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classes on compact manifolds DJl of even dimension 2n. In the construction, we find 
it convenient to use conformally covariant elliptic operators described below. In 
the discussion, we restrict ourselves to even-dimensional manifolds, although the 
corresponding operators can be defined in odd dimensions as well. 

Let OK be a manifold of even dimension 2n, n > 2, and go a Riemannian metric 
on 2Tt. Then, there exists on 971 an elliptic operator P — Pg^ of order 2n, called the 
dimension- critical GJMS operator, constructed by Graham- Jenne-Mason-Sparling 
in |GJMS| . with the following properties: 

P = A" + lower order terms; 

in fact, P has a polynomial expression in (Levi-Civita connection) V and (scalar 
curvature) R, with coefficients that are rational in dimension 2n; P is formally self- 
adjoint ([GZl IFG| ): under a conformal change of metric gi = e^^go, the operator 
P changes as Pg^ — e^^^'^Pgg. 

Given these properties, we can imitate the construction of the GFF in dimension 
2: for /i,/2 G C°°(9Jl), the inner product is defined by 

(/l,/2)p,, ^ / fl{x){Pgj2){x)dv0\g,{x). 

Jm 

Then, this inner product is also conformally invariant. When the metric is 
changed conformally to gi — e^'^go, the volume element changes as 

dvolg, =e2""'dvolgo, 

while P changes as Pg^ = e~^""Pgo. Again, we get the relation 

, /2)Pgj = , /2)PgQ , 

just like in dimension two. 

When n — 2, the dimension-critical GJMS operator 

P4 = + S[i2/3)Rg,go - 2RiCg,]d 

is also called the Paneitz operator. If is flat, then the Paneitz operator is equal 
to A^, hence in M.'^ it is natural to work with A^. However, since Mf^ is not compact, 
we need to consider the operator on a compact domain, in which case we have to 
choose proper boundary operators in order to preserve the conformal covariance 
property. This will be further explored in future work. 

On the compact 2n-dimensional manifold 9Tt, if we construct a Gaussian random 
field using the dimension-critical GJMS operator P, then the covariance function of 
the field is given by the Green's function Gp{x, y) of the operator P. Let d {x, y) be 
the Riemannian distance between x and y on 971. Then, it is known ( |CYllNdllPon| l 
that as d [x, y) 4, 0, Gp{x, y) is asymptotic to — C„ \ogd{x, y) where C„ > depends 
only on the dimension. This is similar to the well-known behavior of the Green's 
function of the Laplace-Beltrami operator A in dimension two. This will become an 
important ingredient in the construction of the random measure on the manifold, 
which we intend to explore in a future paper. 

7. Appendix 

This section contains all the computations with the Bessel functions. We start 
with the Fourier transforms of erf and dcrf , and list all the integral expressions for 
the covariance function of the family {l (/laj) ,1- {hdaf) ■ x £R*,e> O}. 
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Lemma 11. Recall from V2.1]) and S2. 2\) that the Fourier transforms of and daf 
are, respectively, 



and daf {0 = -j/f (0 - -2e~V2 (e |e|) e*("'«)«^ 
Therefore, both and dcr^ are in H^^ (R"') . In fact, for ei, 62 > 0, 



(7.1) [l (/l^xj I (/l^yj 



1 



(7.2) \l [ha^^ j I [hdal 

and 

(7.3) E^ [Z (/id^yj Z 



-1 



r Ji (eir) Ji (e2T) dr, 



27r2eie2 Jo (1 + 



rJi (eir) J2 (e2r) dr, 



27r2eie2 Jo (1 + t2) 



r J2 (eit) J2 (£2^) dr. 



Furthermore, for x,y E M^, x ^ y, and ei, €2 > 0, 



(7.4) 



(7.5) 



and 



(7.6) 



7r2eie2 |x- y| Jo (l + r^) 



r Ji (eir) Ji (e2T) Ji (Ix - y\ r) dr, 



E 



7r2eie2 |x - y| Jo (1 + t2) 



r Ji (eir) J2 (e2r) Ji (|x - y\ r) dr, 



E 



I ( /irf^x ) I ( h^„y 
1 



7r2eie2|x-j/| Jo (l + r2) 



r J2 (cit) J2 (e2T) Ji {\x - ?/| r) dr. 



Proof. Everything follows from straightforward computations in spherical coordi- 
nates in and applications of the following integral expression of the Bessel func- 
tions ( [Wat] . §3.3 ): for every fc > 1 and r > 0, 



Jk{r) 



r(fc + i)r(i) io 

In addition, as we have indicated in Section 2, the asymptotic expansion ([Wat], 
§7.1) of Jfe says that Jk (r) = O (r~^/2) as r — 00. which is sufficient to guarantee 
the convergence of each integral involved in (|7.1l) - (|7.6p 



□ 



It will be convenient to recognize that all the covariance functions involved in 
the previous lemma, i.e.. 



E 



w 



I { h„x ] I \ h 



E 



w 



and E 



w 



{hda-^l{ha„y^ 
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are continuous in all variables x, y G M'* and ei, £2 > 0. In fact, 



I K. ]I h 



d€2 



and 



E 



w 



Tih. 



-E 



I {ha=.^ ^ I {h„y^ ) 



deide2 

These simply follow from the dominated convergence theorem and the fact that 
both Jk {r) and (r) /r are bounded on r € (0, 00) for every A: > 1 . 

Proof of Lemma \^ The proof of all the formulas (I2.3p to ([23]) is based on the 
following integral formulas of Bessel functions which can be found in |Wat| . §13.53, 
pp 429-430: if a > 6 > and p > 0, then 



(7.7) 



r2 p2 



Ji [ar) Ji (6t) dr — Ki (ap) Ii (bp) : 



if a > 6 + c and p > 0, then 
1 



(7.8) 



Ji (ar) Ji (6t) Ji (cr) dr = p-^/fi (ap) h (bp) h (cp) 







+ p^ 

We hereby provide an alternative proof of these two formulas, and complete the 
computations in Lemma [T] For a fixed p > 0, we define the function 



B{a,b) 



1 



- Ji (ar) Ji (br) dr for a,b > 0. 



2Tr^ab Jq + p'^ 

Given a > 0, we observe that if 6x is the point mass at x G M'*, then the following 
integral is finite: 

4 



+ a Id 



27r2a|a;| 



r Ji (ar) Ji (|x| r) 



dr — B (a, |x| 



In other words, 



Therefore, in the 



p^ + r" 

But this integral is also "formally" equal to ^(p^ ^ ^) ^ "'ai '^a 

((p2-A)-V^) (x)==i3(a,|:z:|) 

is a point-wise defined, radially symmetric function in a; G IE 
sense of tempered distribution, 

(p2-A)i?(a,M)=a°, 

which, when written in spherical coordinates, implies 

p^ - dl - B (a, 6) = for aU < 6 7^ a. 

The above is a Bessel-type ordinary differential equation, all the solutions of which 
are in the form of 

^ Ki {bp) h {bp) 

^ ' b ^ ' — b — ' 

where Ci and C2 are two functions only depending on a. Without loss of generality, 
we can assume b < a. If one examines the behavior of B (a, b) when b is close to 
zero, then one finds that Ci (a) = because bB (a, b) converges to zero while Ki {bp) 
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blows up as 6 J, 0. On the other hand, with 6 > fixed, one can apply exactly the 
same arguments to see that B (a, b) also satisfies 

(^p^ -dl- ^d}j B (a, 6) = for all a>h. 

Hence, C2 (a) must be in the form of 



a a 

for some constant C and C". This time, the boundedness of aB{a,b) as a f 00 
implies C" = 0. 

Thus, the only thing left is to determine the constant C. To this end, we observe 
27r^a^B (a,a) — / — 7r^Ji{u)du 



+ a^p 
^ Jl(u) 







u 



du as a I 0. 



However, one can easily verify that 

d_ ( Jl(u)^Jl{u) \ ^ J?H_ 
du\ 2 ) u 

So limaj^o 27r^a^i3 (a, a) — i. Meanwhile, lima^o A (ap) -^1 (fl-V) = \: which implies 
C = Therefore, 

B (a, b) ^ --\-rIi [hp) Ki (ap) for a > & > 0. 
For the formula (17.81). we define 



1 /■°° 1 

C (a, 6, c) = / ^5 ^ Ji (ar) Ji (&t) Ji (cr) dr for a, 6, c > 0. 

■n'^abc Jq + p"^ 

Assume a > b + c. One can verify, by direct computations inside the integral signs 
and the dominated convergence theorem, that 

p'^C (a,b, c) - — ^C(a, &, c) — C(a, 6,c) = for < c < a - 6, 

oc' c ac 

and lim C (a, b,c) ~ B (a, 6) . 

c4,0 

Similarly as above, one has 

B(a,b)^^— . 

cp TT^abcp 

Thus, (fTT)) and ((7?5| are proved. 

Given (|7.7p . notice that 2 ^^o-j^^ I ('^'^fa) '^^'^ computed by applying 
the operator -^-^\p=i to both sides of (|7.7I) with a = ei V £2 and 6 = ei A £2. 
Then fi:om there, based on the earlier observations, the complete expression for the 



C (a, b, c) = -^±^B (a, b) = ^1 (^P) ^1 (^P) 
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covariance matrix in the concentric case can be obtained by taking derivatives in ei 
and 62 accordingly. The detailed computations are as follows. When ei > £2 > 0, 



47r^ V ^2 ei 



Thus, 



E 



de 



9 



E 



w 



_-l ( K[{e^) i^i(ei)/ei \//i(e2)/e2 (/i (£2) /e2)' 



47r2 V K"{ei) (i^i (ei) /ei)' 



Besides, one realizes that 



Ki (61) V 



ei 



i^2 (ei) ^/l(e2)^'_ /2(e2) 



ei 



£2 



£2 



and the formula (|2.3p follows. 

The non-concentric case is very similar. Given (j7.8p . we assign a = ti in case (2) 
and a — \x — y\u\ case (3). By a similar procedure, i.e., applying ■^-^\p=i to (|7.8p 
and taking derivatives in ei and 62, we will be able to compute the non-concentric 
covariance matrix in either (2) or (3). To be specific, when ei > |a; — j/| -I- £2, 



E 



w 



equals 



1 d 
2p dp 



\p=i C (£1,62, \x-y\) 



-1 



, /i (£2) /i - y|) , gi(£i) ,, , , /i(|a:-2/|) 

- , (£ij ^ \ (£2) — 7— 

27r^ V €2 \x-y\ £1 F-2/1 



^ ( K, (£l)/l(£2) ^, _ |. j^i (£i) /i (£2) h {\x ^ y\) 
2tt^\ £1 £2 ^'^ £i £2 \x~y\ 



27r^ V £2 F-yl £1 F-?/| £1 £2 



-^1 (£2) /£2 

(£2) 
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when \x — y\ > ei + 62, ^ (^^<^fj j ^ ('^'^=2 j ^'^^^Is 
1 d . 



C {\x - y\,ei,e2) 



Ki {\x-y\) 



27r2 \ ei £2 \x - y\ ei £2 

1 (KA\x-y\lh{e^) , , K,{\x-y\) ^,^ h 

^ — i i -'i '.'^2) H j j — h (eij — 

2tt^ V F-y| ei \x-y\ £2 

-1 ^ „/i(£i)/i(£2) , K,{\x-y\) h{t,)^, ^ K^{\x-y\l h 

Zlf MlA /((£,) ^''"''-^^'""''^ nfer V^i(^2)/£2 



The rest is straightforward. Thus we have finished the proof of Lemma [1] □ 

Before continuing, we point out that similar computations can be carried out 
in higher even dimensions R^n with n > 2. In fact, the formulas (|7.7I) and (I7.8P 
remain true if one replaces Ji by Jn-i in the integrands, times a factor of r^"" 
to the integrand of (|7.8p . and replaces i^i by h by 7„_i respectively in the 

results. Therefore, one can use these modified results to compute the covariance 
function of the Gaussian family : x £ M^",£ > 0} (as defined in Section 

6.1) by applying the operator i—^j^j \p=i to the modified version of (|7.7p and 
(|7.8p . The rest follows similarly as above. 

Next, we want to provide details in the deriving the formulas for /i^ (|2.6p and 
G (e) (1^ as weU as the results (P?5| - (P7TU| in Theorem!^ It's an easy matter to 
check that for every £ > 0, 

det B (£) = £-1 (If (£) - Jo (£) I2 (£)) > 0, 

where we applied the Bessel function identities ( |Wat) . §3.71) /( (£) = — + ^0 (e) 
and /(' (£) = + h (e) . Therefore, 

B-i (,) _ 1 f (e) - ^2 (e) h (e) - elo (e) 



/?(£)- Jo (£) /2 (£) V -^2(e) hie) 
Recall that = B^^ (£) Vf, when computed explicitly, 

jj. ^ 1 f {eh (e) - h + ih (£) - £/o {e))l{hd..) 

'"^ lUe)- lo (e) I2 (e) ^ ~h (e) X (/^.j ) + h (e) X {hda^ ) 

and if C = (I7 l)^i then fi^ — fi {e) erf + /2 (£) da^ where 
(7.9) A (£) = -27- J , .J , . and A (£) 



H (£) - /o (e) hie) H (£) - h ie) h ie) ' 

from which one sees that /if has the "right" limit as £ 0. In addition, we can apply 
more Bessel function identities: 

1[ ie) = \h ie) + h ie) , h ie) - h (e) = \h (e) , 
K[ (£) - ii^i (£) - K2 ie) and h (e) i^2 (e) + ^2 (e) i^i (e) = ^, 
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to write down B ^ (e) A (e) explicitly as [if (e) — /o (e) h (e)) ^ times 



-(1 + e-) 
h{e)Ki{e)+l2{€)Ko{e) + e-^ 

Therefore we have: 

(1), given X gM.^ and ei > e2 > 0, 



/i(e) K, (e)+/2(e)ifo(e) + e- 



E 



x(/z^yjx(/i^yj] =G(6)^ 



1 

1 



C^B-i(e)A(e)C(with C^ = (l,l)) 
2Ji(e)ifi(e) + 2J2(e)i^o(e)-l 



7f(e)-7o(e)/2(e) 



(2), given a;, y e M^, x ^ y, and ei, e2 > with ei > |a; — t/|+e2, E^ 
equals ^ (if — /0/2) ^ (ei) times 
C^B-i (ei)A(ei)C(|x-2/|)C. 



= ( {hKi + hKo) (ei) - 1 (/li^i + I2K0) (ei) ) 
'2/i(|a:-y|) 



\x-v\ 



+ -f2(|x-2/|) (/iXi + /2ifo) (ei) 



+ \l2{\x-y\) 



=Io {\x - y\) i^ihK, + I2K0) (ei) - - j + -h {\x - y\) ; 
(3), given x,y G W^, x j^y, and ei,e2 > with \x - y\ > ei + e^, 



E 



w 



1 



1 
2"^ 
1 

271"^ J \ \x — y\ 
Ko{\x-y\). 



f2Kii\x-y\) 



-K2i\x-y\) 
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